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I ARTH S GRAVITY I Ihl.l) MAPPIN(i 


RFOUIRI Ml NTS AND CONCI PT 
P. (). Vonbun uml W. I). Kahn 

ABSTRACT 

Tlic objective of this paper is to discuss a future sensor (gradiometer) for mapping the earth’s 
gravity field to meet future scientific and practical requirements for earth and ixeanic dynamics. 
These are approximately tO.I to 10 mgid over a block size of about 50 km and over land and an 
ocean geoid to I to 2 cm over a distance of about 50 km. 

To achieve these values requires a gravity gradiometer with a sensitivity of approximately 10"^ 
I'll in a circular polar orbiting spacecraft with an orbital altitude ranging 160 km to 1 80 km. 
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Rl OUIRI MI NTS AND CONO PT 

INTRODUCTION. (iRAVITY ANOMALY RI OUIRI Ml NTS 

TIte observation of spacecraft orbits started a revolution in the determination of the “shape" 
of our planet. ( Ref. I ). As time has progressed, so have the demands for gravity information. Spac>*- 
craft orbital perturbation an Ivses were and still are the principal means for deriving the earth gravity 
field, llie best gravity field model published in the open literature to date may be the I OH 
(Ref. 2), a field with a block size resolution of about SSO km (S* x S°) and an error of about 10 
mgal on the average. 

It is well known that the earth's gravity field is a result of the distribution of its mass and thus 
is the fundamental base needed to construct an accurate geoid. Tlie latter in turn is the ..leady state 
surface of the ocean if there were no tides, currents, wind pile-ups and atmospheric pressure changes. 
On the other hand, it is these changes which can tell us about the ocean currents if a reference surface 
(geoid) is known to an accuracy of I to 2 cm and a spatial resolution of approximately 50 km and 
larger (Ref. 3). It is this small resolution which dictates a new approach to the detennination of the 
gravity field. 

At present the best approach to increasing our knowledge of the earth's gravity field is the use 
of Satellite-tivSatellite Tracking (SST), which was started at Goddard Space Flight Center (GSFC) 
during the late 1960's. These techniques demonstrated experimentally with ATS-6 and Nimbus E, 
with ATS-6 and GFOS-3 (Ref. 4), and finally with ATS-6 and Apollo-Soyuz (Ref. 5) which can be 
considered the first demonstration of a low-flying type gravity satellite mission. Tliis effort has now 
culminated in GRAVSAT-A (Ref. 3 and 6), a possible 1983 new program start candidate as an Appli- 
cations Spacecraft essential for future studies of the dynamics of the solid earth and the oceans. 
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I1u- scientific re(|iiirenieiits tor tiilure stiiilick of earth and iKcan dynamio can be siimmarued 
as follows ( Ref. 3): 

• Solid l-arth Physics: tO. I to I.O iniial; 50 to bO kin resolution 

• Ocean ('irculation; t I to 2 cm; 50 to bO km resolution 
riiese requirements can he compared with our present knowledite; 

• Solid l-arth Physics: lO mgal; 5(K) km resolution. 

• Ocean Circulation: 2(W cm gcoid; 5(K) km resolution. 

• loKal Ocean areas: 20 mital; 1 00 km resolution. 

50 cm; l(X) km resolution. 

MISSION OHJI (TIVI S 

Pie objectives of the (.iRAVSA f-A mission to satisfy the above requirements can now be stated 
as follows: 

• l*rovide an accurate gravity field to studv the structure of the crust and upper mantle of the 
earth. 

• Provide an accurate global geoid to be used as the referem.c surface for studies of ocean 
currents. 

CON( I Rf 
Sensor Systems 

As has been mentioned, the use of SST will jvnnit resolution of gravity field anomalies with an 
accuracy of about ± I mgal for each I“ X I" area of the globe, the present goal and capability for 
CiRAVSAT-A. However, this seems to be the limit of the SST technique as shown below. 

If anomaly resolution is to be universal for a block size of approximately 50 X 50 km another 
“sensor” must be developed and used. A gravity gradiometer with a sensitivity of 10~^ l:U can 
achieve these more stringent requirements. 



In the I'ollowinit. the relationship between the two possible "sensors," tor mapping the gravity 
field, i.e„ SST and a gradinmeter ( Kefs, 7. 8) is shown. The reason for considering the gradiometer 
is that it can resolve the sl.ort wavelength components of the earth's gravity field up to degree and 
order 4(H). 

In order to estimate and compare the expected "Signals," we shall analyze their respective power 
spectra ( Kefs. 7, H, ‘J). The earth’s gravity potential is given by (e.g. Ref. ‘M 

Pj^(sin0)|r,j^cosmX + sininM. (I) 

L 8" ’ tn-O 

where 

r IS the geocentric s:itellite distance, 

0 is the sub-satellite geocentric latitude, 

X is the sub-satellite east longitude, 
a is the mean radius of the eartl., 

I\^(sin0) aie the fully normalized as.s(Kiated Legendre polynomials of degree B and order m, and 
Cjn, . S^^ are the fully normalized gravitational cwfficients. 

These- power spectra will nvvw be derived for the SSI and the gradiometer systems. 

Satellite-to-Satellite Tracking System (SST) 

The SST system measures the range rate between two spacecraft (Ref. 6) resulting from the 
difference of gravity acting upon the two spacecraft. The total energy of a spacecraft in orbit is 
given by; 

L * Vi mv* + T (*) 

where v is its speed of the spacecraft, m its mass and T is the disturbing potential of the earth. The 
variation of T is then 

6T = -v6v. (3) 
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i‘foni equation ( I ) withtMit the cential force tenn, together with equation (3i we obtain 


6v ■ V 



^(sin0)(Cj^cosmX + S^^sinmX). 


(4 


*•2 m »0 

The veltH'ity power spectrum of 6v is now obtained by squaring equation (4) and integrating over 
the unit sphere: 




(5 


ni*0 


where Kaula's rule (Kef. 10), the summation term can be replaced with 




Ilie S(|uare riK)t of the power spectrum, that is the expected velocity signal for each wave number 
reads: 


where 


(ST trj 


c« + 



(( 


h » spacecraft orbital height, 

a K mean radius of the earth (i.e. nearly circular orbits are assumed for the missions 
discussed), and 

(>M B gravitational parameter of the earth. 

introducing (iM ■ 4 X 10*® cnr^/s‘, assuming K > 40, and convertinr. to /im/s results in: 

S, (6v) » :.82 X 10’ (a + h)"'^ ( ^ (' 

for a and h expressed in cm. 

These spectral values are described in more detail and compared with experimental data in the 
next section. 
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(OMPARISON WITH THI APOLLO SOYUZ HXPl RIMLNT 


Salollite-to-SatcIlitc Kanjtc Kjlc SignaU 

As was itiued earlier, an actual •alell'te-ttvsalellite tracking (SSI'I experiment ( Kef. 5) has 
demonstrated that the recovcr>' of gravity anomalies is feasible. In the later Apollo experiment the 
measurement noise of the ATS-b to Apollo link (higli-low concept) was determined to be about 300 
pm/s. Tlie next step it to relate this experimental value of 300 ^m/t to the value predicted by equa- 
tion (6). Tliat is. for the Apollo orbital altitude of 230 km and a 550-km horizontal resolution 
(C ■ 3b), the signal amounts to 

Tins, however, represents a global value. Since the experiment was conducted over a hKal region, a 
gravity signal increase was to be expected. 

('•ravity Anomaly Uncertainties 

One of the important quantities to be estimated for SST techniques is the uncertainty attribut- 
able to the technique as was discussed before. A theoretical estimate of the errors in terms of the 
nrcovered gravity anomalies can only be obtained by simulation studies which are in general tedious 
and costly. However, if experimental values are available, simple scaling of these values can be applied 
for error estimation. In this case, the Apollo- Soyuz errors will be used. In order to obtain the scaling 
law, we must express the gravity anomalies in terms of the disturbing potential T (Ref. *)). that is 

3T : 

5g * - T T. (9) 

3r r 
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In spherical harmonics, equation (<)) reads as I'ollows: 


CM *""* 

—7 "(^)* '(Cj^cosinX •»• Sj^sinmX). 


( 10 ) 


««2 m »0 


An estimation of this increase (noise decrease) is obtained by adjusting the global value, i.e,, by multi* 
plying the signal by the ratio of the global surface area to the local experiment area. Tliis problem is 
treated in detail by I-. II. (iaposchkin ( Ref. 12). Tor the Apollo experiment this value is about 3 
(Ref. 5). 

Applying this priKess to equation (8) yields a signal of about 320 ^m/s. As.suming further that 
there were an average two cross-overs per 5X5° block and that the signal is noise limited, the signal 
is a factor of V 2 larger, that is 450 ^m. It is assumed that the signal used in the computation is a 
composite, in the least-square sense, of all the signals detected over the experimental area. Since the 
predicted single spectral signal is about 140 /am s and the “useful” signal was about 450 /am, the com- 
puted signal-tivnoise ratio was found to be about 3, making a gravity estimation possible. Tlie anal- 
ogous expression to equations (4) and (5) for the power spectral value of 5g is 




2c c 
m ■ 0 


( 11 ) 


Simplifying equation ( 1 1 ) by assuming C » I and replacing the sum of the normali/.ed ax'ITicients 
of the harmonic expression by Kaula's lule (Ref. 10) one obtains the gravity anomaly signal 


S ( 6 g) » l—]x lO-!' cm/s 2 . 

(a + l.)^C^Va + h; 


( 12 ) 


Using equation (I 2), one can now form the needed ratios representing the desired scaling law': 


SjjiSgi) a + hj y(J,y ^a + h,j \a + h,j 


(13) 
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when.* tuhschptk I and 2 rrftrr to Apollo and (travul rokpcctivcly, 
h, - 230 X 10^ cm. 
h, ■ 160 X 10* cm, 

5j *36. corrcspondintt to a (550 X 550) km ^ block «/e. 
t, ■ 200, corrcspondinit to a ( 1 00 X 100) km^ block size, 

Vj| I6K| ) ■ 10 miial (R.M.S. associated with recovered (550)* km* anomalies from Apc)llo 
based on companson with (il-M lOB, 

I 6k-») • predicted rms for (Iravsat ( 100 X |(X))km* anomaly, 

N| * number of orbit revolutions available for Apollo derived anomalier <Nj ■ 2), 
N, * number of orbit revolutions expected for (iravsat derived anomalies (Nj ■ 4h), 
a ■ mean earth radius ;6.4 X 10* cm) and 
CiM ■ 4 X 10*® cm*/scc*. 

We now form the ratio of et|uation (6) to Equation ( I 2): 


I0(a -t-h)*/^ 

C’ombinint! (13) and (14) then yields 

I0(a + h,)*/2/^ + 


/itn /s/mi:al 


(14) 


-K, 


' 2 , N . n '/2 


S,j(6vj) 


(l^) U) 


(CiND'/^g^ V2/ 

Using the values specified above in equation (15), one obtains an estimate of the precision re- 
quired for the (iravsat-A sensor, that is 


Sg (6v,) * 0.8 /im/s. 


This value agrees quite well with the precision specified for Ciravsat-A. I ^m/s ( Ref. 3). 

From equation ( 15) it can be seen that only a moder’’te improvement or reduction of error m 
6g can be obtained by lowering the (iravsat altitude hj (Nj is increased correspondingly). For 
example, Ic vering the orbital altitude from 160 km to 150 km changes the required sensitivity by 
only a factor of 1.4; but such a step introduces serious engineering problems (more fuel, tighter 
control, thermal problems, etc.). 
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(iRADlOMI Tl R SYSTI M 


llsinft iin appitMch similar to that used for salcllilc-to-satcllitc tracking (SST) we can obtain 
spectrum for the gravity gradient wnsed by a super cmiled gradiometer system as follows; the 
second derivative of the disturbing potential T with respect to the radius. d^T/Ar^. is the radial com' 
ponent of the gravity field tensor. This com|H)iient is the largest in magnitude of the nine compi>- 
nents of the gravity tensor. Torming the ratio of the radial to the along-track and crosstrack 
components of the field tensor one obtains (Rtf. 7). 


V'J (grad) V* igrad) 

*rid _ *ti<J 


V* (grad) V* (grad) 

cruM <k>ns 


i 4. 


(17a 


Tlie other ratios, involving the cross-components are very much smaller. I'luis the spectrum of the 
main component is used here for estimating the expected signal for the gravity gradiometer. 

Tlie radial component of the gravity field tensor is 
N 


a 2 - 


*T (iM v~< /aV „ f- 

-7 • — ^ ^ (V+l)(t + :)(-j Pj„,(siiv) j(\„,cosmX + S,„,sinmX 


v2 m*ii 

Tlie power s(vctrum by degree K is given by the stpiarc of the sum over all the m’s. (hat is 


';tgrad).^yj (K * I )*(« + 2)- ) . 


m • 0 


l or S >40. the signal in simplified form is: 


(17b 


(IK 


Sj,(grad) ■ (2 X 10*“’)'' 


C.M 

(a + h)^ 



(t)l/2 


cm ;s* 
cm 


(IKa 


Introducing the value for (iM. we obtain the “signal" by degree in terms of I II; 


(ll-U ■ lO"** cm/s^/cm). 


Sj(grad) 


5.b4 X 10^* 
(a + h)^ 



(C)'' (KU). 


(19 


For the proper numerical valuet of the expected iianalx one ohtaint from (7) and (IKa) 


reiiiectively 


and 


Sj8v)A(II X («)”*'* Mm/t 


S(grad)*C X lO*)" 


c-^x-^r 


(«)'* I II 


COal 


COb) 


where 


a ■ the mean radius of the earth, 4X10^ km 
h ■ the orbital heiglit of the spacecraft, IbO km and IKO km, 
g • acceleration of gravity, d81 cm/s^, and 

V ■ the degree of the spherical harmonic expansion of the earth gravity field. 

As can be u'en, the gravity fi- d attenuation factor (a/a r- h)^ it the same in each cate. Hie 
quantity (a/a-rh)' for h < a, as is always the case for gravity type missions, is neat unity (0.^)5). What 
is interesting, however, is that the resolution C (degree of the harmonic expansion) has a drastically 
different effect on the expected "gravity signals'* in favor of the gradiometer for (i > ISO. In terms 
of block si/e b: where 

^ 2ira . :0.0()() , 

b ■ km , 

2K « 

I'igures I and 2 depict the expected gravity spectral signals as sensed by the SSI and gravity 
gradient sensors for orbital altitudes of l(>0 and IKO km res|vectively. A relationship can be obtained 
*'rom fonning the sipiare rmit of the ratio of equations ( 5 1 and ( I K) that is 

j ’ (-^) (,p) S.,(erad) ( 


Assuming h| * lu and equation (21 ) reads 

f(a + h)*l \\ 


S((6v) cm/s ' 


S (grad) (cm/s^/cm) 
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Iniroiluciny fttnii ami l-U. <mc obtaini 


S,(6v)^m/» a S^(grad) 



I U 


(21a) 


I'or h ■ |60 kin ami 8 ■ 200, 


Sj(6v)pin/«* 137 (grad) HI. (2lh) 

One »iin/t lhu» corrciiponds to about 7.10"^ Ml. 

In order to make a meaningl'ul a»M*Mment of the above spectral signals the wave number spectrum 
it integrated for different intenals of degree C. 

Dm provides a measure of the nuHlelled signal alx>ve a certain degree. Integrating i (|uat;ont 
20a and 20b yields 


Si ■ r * Sjifiv) di 

A, 

l-or the SST concept and 

S j^(grad)" I Sy(grad)dC 

* I 

for the gradiometer apparent. 

Tables I & 2 show the signal .drength for selected intervals of degiee C. 


Table I 

SSI Wave Number Spectrum 


h (km) 

«1 


Si ,*.j(fi' )/im/sec 

160 

36 

100 

3.6 X 10-^ 


100 

180 

2.3 X 10- 


180 

360 

16.4 X lOO 


360 

1000 

7..SX 10-2 

180 

36 

100 

3.1 X 10-^ 


100 


1.6 X 102 


180 

360 

8.6 X lOO 


360 

1000 

2.3 X 10-2 


(22a) 


(22b) 













TiiNc 2 

(iriiilioniclcr Wave Number Spectrum 


h (km) 


* 

S^li^lgrad) hU 

160 

36 

100 

2.05 


100 

180 

0.67 


IKO 

360 

0.14 


360 

1000 

2X lO-J 

180 

36 

100 

1.68 


100 

IKO 

0 45 


IKO 

360 

007 


360 

1000 

6X 10-^ 


Tl«e vaiue« lahulatcJ show that comparatively large signals exist compared to the systems men- 
tioned ahtwe (i.c.. I ninlt, 10"^ I'Ul. 

l or instance lor I KO ^ C < 3(>0 the signai is about 16 pm s and I4(X) X 10“^ HJ for orbital 
height of 160 km. 

11m demonstrates that the high frequency components of the gravity field exhibit a signal to 
noise ratio of 14 to I and 1400 to I res|>ectively. Tins underscores the power of a gradiometer type 
space mission. 

A.S can be seen from Tigures I and 2. and Tables I and 2. a gradiometer with a sensitivity of 
10"^ ITI can achieve a resolution of around 50 to 60 kin. 

Such a system is in the development stage under a NASA grant at the University of Maryland 
(Figure 3). Figures 4 and 5 are extensions of F'lgures I and 2, respectively. These show clearly the 
resolution limits of both spaceborne systems even if the systems* sensitivities are improved by several 
orders of magnitude: the last number representing the ultimate limit for a supercooled (4*K) gradiom- 
eler sensor according to Dr. Paik ( Ref. 1 1 ). 

As can be seen a sensitivity chan(% of the gradiometer of four orden of magnitude ( 10“^ to 
10“^ FU) results only in a resolution change of a factor of two (from 60 km to 30 km). This in 


II 











cncncc means that a gruJiomelcr with a sensitivity of, say. I0~^ TU constitutes a practical limit. 
Any lurtlwr improvement in the sensitivity will prtxiiice neidigibir results as far as spatial resolution 
is concerned. 

We next discuss the effect of the gradiometer precision on the gravity anomaly error. We first 
form the ratio of eipiations ( I Ml and (III; that is 


- - • I cm . 

.S^(grad) \ t / 

Then, scaling to IT' and ingal yields 

SjiAgi (ingal) * 10"* ^^^-^jsyigrad lITI, 
where V! is nondimensional and a and h are in cm. 


( 22 ) 


(23) 


I'quation ( 23) states clearly that an increase in the gradiometer precision diws suhslantially de> 
crease the error in the estimated gravity field, l or example, using gradioiiKter orhithig at l(>U km, 
and a I* X I* gravity anomaly (K » 2(X),t, one obtains from (23) 

S^(6g) mgal « 3Sj(grad) IT' (24) 

A gradiomc*er with a precision of 10“^ yields a gravity anomaly enor of 3 x 10“^ mgal. lliis 
would mean that using such an instrument one cun resolve gravity field anomalies of about 50 X 50 
km to less than one m(uI. Practically, this value will probably be less than one mgal since computa- 
tioi'.u difficulties (high correlation etc.) will certainly arise. Only a successful but computationally 
costly simulation will resolve this question. | 

('ONCLL'SIO:' 

In summary it can be stated that » further improvement in sensitivity and resolution of the 
earth gravity field can be realised using a supercooled gravity gradiometer instead of satellite-to- 


satellite tracking system. 



The former sentor provides in situ measurements where the latter pro ides derived field param- 
eters. A spacchurne polar orbitinn gradiometer with a precision of I0~^ hU operating at an altitude 
of IhO km can map the gravity field to a precision of say less than one mgal and a horizontal resolu- 
tion of about 50 km. This further corresponds to a geoid height error to the suheentimeter level, 
wtiich in turn satisfies almost all future ocean dynamics requirements. 
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